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The nucleation process of quark matter in cold (T = 0) stellar matter is investigated using the 
microscopic Brueckner-Hartree-Fock approach to describe the hadronic phase, and the MIT bag 
model, the Nambu-Jona-Lasinio, and the Chromo Dielectric models to describe the deconfined 
phase of quark matter. The consequences of the nucleation process for neutron star physics are 
outlined. Hyperonic stars are metastable only for some of the quark matter equations of state 
considered. The effect of an hyperonic three body force on the metastability of compact stars is 
estimated, and it is shown that, except for the Nambu-Jona-Lasinio model and the MIT bag model 
with a large bag pressure, the other models predict the formation of hybrid stars with a maximum 
mass not larger than ~ 1.62 Mq. 
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I. INTRODUCTION 

Quark matter (QM) nucleation in neutron stars has 
been studied by many authors both at zero [H-Q and 
finite temperature [ll| - |l6| . due to its potential connec- 
tion with explosive astrophysical events such as super- 
novae and gamma ray bursts. In all these works, the 
hadronic phase was described using phenomenological 
models, such as e.g., the well-known relativistic mean 
field (RMF) model based on effective Lagrangian den- 
sities where the baryon-baryon interaction is described 
in terms of meson exchanges 17). Among the different 



RMF models, one of the most popular parametrizations 
is the one of Glendenning and Moszkowski [l8| of the 
non-linear Walecka model which have been widely used 
to study the effect of the hadronic equation of state (EoS) 
on the QM nucleation process. In particular, the effect 
of different hyperon couplings on the critical mass [l], Q 
for pure hadronic stars (HSs, i.e. neutron stars in which 
no fraction of QM is present) and the stellar conversion 
energy [l9| was studied in Ref. Q- It was found that 
increasing the value of the hyperon coupling constants, 
increases the stellar metastability threshold mass and the 
value of the critical mass, thus making the formation of 
quark stars (QSs, i.e. hybrid stars or strange stars de- 
pending on the details of the EoS for quark matter used 
to model the phase transition) less likely. In that work, 
the hadronic phase was also described using the quark- 
meson-coupling model concluding, in that case, that 
the formation of a quark star was only possible with using 
a small value of the bag pressure. In all these works the 
MIT bag model [2l[ was used to describe the quark mat- 
ter phase. In a recent work [l(| , two models that contain 



explicitly the chiral symmetry were applied to describe 
the quark phase, namely the Nambu-Jona-Lasinio (NJL) 
model [13] (see also [23l. l2~j] ) and the Chromo Dielectric 
model (CDM) [25l.l26l|. It was shown there that it is very 
difficult to populate the quark star branch with the NJL 
model and, therefore, all compact stars would be pure 
hadronic stars in that case. On the contrary, with the 
CDM, both hadronic and quark star configurations can 
be formed. 



In the present work we study the nucleation of quark 
matter using an hadronic EoS based on microscopic cal- 
culations. In particular, wc employ two hadronic EoS 
based on microscopic Brueckner-Hartree-Fock (BHF) 
calculations of hypcrnuclear matter. The first one (here- 
after called Model 1) is the recent parametrization pro- 
vided by Schulzc and Riiken |27l which uses the Argonnc 
V18 nucleon-nucleon (NN) [28| supplemented by the mi- 
croscopic three-body force (TBF) of Ref. [29| between 
nucleons (NNN), and the recent Nijmegen extended soft- 
core ESC08b hyperon-nucleon (YN) potentials [HO] - The 
second one (hereafter called Model 2) is based on our re- 
cent work of Ref. [3l[ where we used the Argonne V18 
NN force and the Nijmegen soft-core NSC89 YN one [32 1 
in a microscopic BHF calculation of hyperonic matter 
supplemented with additional simple phenomenological 
density-dependent contact terms, that mimic the effect 
of NNN, NNY and NYY TBFs, to establish numerical 
lower and upper limits to the effect of hyperonic TBF 
on the maximum mass of neutron stars. To describe the 
quark phase, in the present work, we use the three differ- 
ent models already mentioned, the MIT bag model 21 
the NJL model [12] and the CDM model [2j|. 



The paper is organized in the following way. In sec- 
tions |n] and Unl wc briefly review the BHF approach and 
the main features of quark matter nucleation in hadronic 
stars, respectively. Our results are presented in Sec. IIVI 
Finally, a summary and the main conclusions of this work 
are given in Sec. [V] 
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II. THE BHF APPROACH 

The BHF approach is the lowest order of the 
Brueckner-Bethe-Goldston (BBG) many-body theory 
[33[ | . In this theory, the ground state energy of nuclear 
matter is evaluated in terms of the so-called hole-line ex- 
pansion, where the perturbative diagrams are grouped 
according to the number of independent hole-lines. The 
expansion is derived by means of the in-medium two- 
body scattering G matrix. The G matrix, which takes 
into account the effect of the Pauli principle on the scat- 
tered particles and the in-medium potential felt by each 
nucleon, has a regular behavior even for short-range re- 
pulsions, and it describes the effective interaction be- 
tween two nucleons in the presence of a surrounding 
medium. In the BHF approach, the energy is given by 
the sum of only two-hole-line diagrams including the ef- 
fect of two-body correlations through the G matrix. It 
has been shown by Song et al. [34| that the contribu- 
tion to the energy from three-hole-line diagrams (which 
account for the effect of three-body correlations) is min- 
imized when the so-called continuous prescription [35j is 
adopted for the in-medium potential, which is a strong 
indication of the convergence of the hole-line expansion. 
The BHF approach has been extended to hyperonic mat- 
ter by several authors [27], HH, [3(| . The interested reader 
is referred to these works for the specific details of the 
BHF calculation of hyperonic matter, and to Ref. [33| for 
an extensive review of the BBG many-body theory. 



III. QUARK MATTER NUCLEATION IN 
HADRONIC STARS 



The conditions of phase equilibrium, in the case of a 
first-order phase transition [37j, are given by the Gibbs' 
phase rule, which in the case of cold (T = 0) matter can 
be written as: 



thus v. 
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This time is many orders of mag- 



P H = P Q = P a 
where 



Vh(Po) = Vq(Po) (1) 



m = MQ = — 2 

n H n Q 

are the Gibbs energies per baryon for the hadron (H) 
and quark (Q) phases, respectively, and the quantities 
e H (cq)j Ph(Pq), an d uh&q) denote respectively the to- 
tal {i.e., including leptonic contributions) energy density, 
total pressure, and baryon number density of the two 
phases. Above the transition pressure Pq the hadronic 
phase is metastable, and the stable quark phase will ap- 
pear as a result of a nucleation process. Quantum fluctu- 
ations will form virtual drops of quark matter. The char- 
acteristic oscillation time Vq 1 of these drops, in the po- 
tential energy barrier separating the metastable hadronic 
phase and quark phase, is set by strong interactions, 
which are responsible of the dcconfmcmcnt transition, 



nitude smaller than the weak interaction characteristic 
time (T W eak ~ 10~ 8 s), consequently quark flavor must 
be conserved forming a virtual drop of quark matter. We 
call Q* -phase this deconfined quark matter, in which the 
flavor content is equal to that of the /3-stable hadronic 
phase at the same pressure and temperature. Soon after 
a critical size drop of quark matter is formed, the weak 
interactions will have enough time to act, changing the 
quark flavor fraction of the deconfined droplet to lower its 
energy, and a droplet of /3-stablc quark matter is formed, 
(hereafter the Q-phase). 

This first seed of quark matter will trigger the con- 
version [r| [H, HI] of the pure hadronic star to a quark 
star. Thus, pure hadronic stars with values of the cen- 
tral pressure P c higher than Pq and corresponding masses 
M > Althr = M(Po), are metastable to the decay (con- 
version) to quark stars [H-Q. The mean lifetime of the 
metastable stellar configuration is related to the time 
needed to nucleate the first drop of quark matter in the 
stellar center, and it depends dramatically on the value 
of the stellar central pressure. 

As in Refs. [H-i], we define as the critical mass M cr 
of the hadronic star sequence, the value of the stellar 
gravitational mass for which the nucleation time of a Q*- 
matter droplet is equal to one year: M cr = Mhs{ t — 
1 yr). Pure hadronic stars with M > M cr are thus very 
unlikely to be observed. M cr plays the role of an effective 
maximum mass Q for the hadronic branch of compact 
stars. 

In a cold and neutrino-free hadronic star the forma- 
tion of the first drop of quark matter could take place 
solely via a quantum nucleation process. The basic quan- 
tity needed to calculate the nucleation time is the en- 
ergy barrier separating the Q*-phase from the metastable 
hadronic phase. This energy barrier, which represents 
the difference in the free energy of the system with and 
without a Q*-matter droplet, can be written as pit |40T] 

U(K) = ^nn Q , (jig. - m )K 3 + AttuTZ 2 , (3) 

where 71 is the radius of the droplet (supposed to be 
spherical), and a is the surface tension for the sur- 
face separating the hadron from the Q*-phase. The 
energy barrier has a maximum at the critical radius 

7l c = 2o-/[n Q *((i H -Mq*)]- 

The quantum nucleation time T q can be straightfor- 
wardly evaluated within a semi-classical approach 0- 
0, H3| and it can be expressed as 

r q = [v m) N c )- 1 , (4) 

where po is the probability of tunneling the energy barrier 
U (71) in its ground state, vq is the oscillation frequency 
of a virtual drop of the Q*-phase in the potential well, 
and N c ~ 10 48 is the number of nucleation centers ex- 
pected in the innermost part (r < R nuc ~ 100 m) of the 
hadronic star, where the pressure and temperature (in 
finite T case) can be considered constant and equal to 
their central values. 
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IV. RESULTS AND DISCUSSION 

We will now discuss the results obtained with the two 
microscospic hadronic EoS considered and, in particular, 
we will comment whether the possible discussed scenarios 
are compatible with the recent measurement (4lj of the 
mass of the pulsar PSR J1614-2230 with a mass M = 
(1.97 ± 0.04) Mq. 

In Fig. [1] we plot the Gibbs energy per baryon as a 
function of pressure using the microscopic approach of 
Ref. [27j (Model 1) for the hadronic phase and one of the 
following models for the Q*-phase: MIT bag (top pan- 
els), CDM (left bottom panel) and NJL (right bottom 
panel) model. It is interesting to note that the forma- 
tion of the Q* -phase is possible only in the case of the 
MIT bag model EoS with a low value of the bag constant 
(B = 85 MeV fm~ 3 ). In all the other cases considered 
in Fig. [TJ the curve for Gibbs energy per baryon for the 
Q*-phase never crosses the one for the hadronic phase, 
consequently, the hadronic phase will always remain sta- 
ble with respect to the formation of Q*-phase droplets. 
For these three QM models, this result implies that the 
pure hadronic stars (hyperonic stars) described by Model 
1 are stable up to their maximum mass configuration 
M™, = 1.37M . 

By numerical integration of the Tolman- 
Oppenhcirmer-Volkov equations [42j , we have calculated 
the structural properties for pure hadronic and quark 
star sequences. The main results, in the case of Model 
1 EoS for the hadronic phase, are summarized in Table 
HI where we report the maximum gravitational mass 
Mmax (third column), the gravitational threshold mass 
Mthr = M(Pq) for metastable configurations (fourth col- 
umn), and the gravitational (baryonic) critical mass M cr 
(M c b r ) (fifth (sixth) column) for the pure hadronic star 
sequence. Mfi„ (seventh column) is the gravitational 
mass of the hybrid star formed by the stellar conversion 
process of the HS with M = M cr and assuming baryon 
number conservation in the process (i.e.. assuming 
Mj in = M^ r ). Finally E conv is the total energy liberated 
in the stellar conversion. It is interesting to note that 
hybrid star configurations can be obtained with the 
MIT bag model with B = 85 MeV fin" 3 , and with 
the CDM. In the latter case, however, the transitory 
non-/3-stable Q*-phase is not energetically achievable 
(see left bottom panel Fig. [1]). Thus in this case the 
hadronic star sequence is stable up to the maximum 
mass configuration (thus we have no entries in Tab. U 
for the quantities Mthr, M cr , M% r , Mf in and E conv ). 
In the case of the MIT bag model with B = 85 MeV 
fm~ 3 the Q*-matter nucleation is possible and one has 
M cr = 1.272 Mq. The conversion of this star will 
produce an hybrid star with Mfi„ = 1.233 Mq. If this 
object is a member of a binary stellar system, eventual 
accretion of matter from the companion will allow it to 
reach a maximum mass of 1.544 Mq. In this case the 
pulsar PSR J1614-2230 will neither be an hyperonic star 
nor an hybrid star. 



We have also artificialy turned off the hyperonic de- 
grees of freedom in Model 1 EoS, and considered pure 
nucleonic stars. In this case we have M^ x = 2.27 Mq, 
and a critical mass M cr ~ 2.2 Mq for all the quark mat- 
ter EoS considered in Tab. HI In all cases, however, the 
critical mass configuration will collapse to a black hole 
(BH entry in seventh column in Tab. [TJ) and thus the 
hybrid star sequence can not be populated. In this case 
the pulsar PSR J1614-2230 would be an hadronic star 
containing only nucleons and leptons. We next dis- 
cuss the results obtained with the hadronic EoS based 
on our recent work of Ref. [3l| (Model 2) where, as we 
said, a microscopic Brueckner-Hartree-Fock approach of 
hyperonic matter based on the Argonne V18 NN and the 
NSC89 NY forces is supplemented with additional sim- 
ple phenomenological density-dependent contact terms 
that mimic the effect of nucleonic and hyperonic three- 
body forces. In particular, we consider three different 
parametrizations of this model corresponding to differ- 
ent values of the incompressibility coefficient, A"^, of 
symmetric nuclear matter at saturation, and the param- 
eter x, which characterizes the strengh of the hyperonic 
three-body forces: Pi (K^ = 236 MeV and x = 1), P 2 
(A'oo = 236 MeV and x = 1/3) and P 3 (K^ = 285 MeV 
and x = 1). The interested reader is referred to [3l| . 
and particularly to Tables 1 and 2 of this reference, for 
details. The results for this model are shown in Table HT1 
and Fig. [2j As in the previous case, we summarize in 
Table [XXI the main stellar properties obtained with this 
model in combination with the MIT bag and the CDM. 
Note that in this case with the NJL model no transi- 
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FIG. 1: (Color online) Gibbs energy per baryon at zero tem- 
perature as a function of the pressure for the hadronic phase 
(solid lines) and the Q*-phase (dashed lines). Model 1 has 
been used to describe the hadronic phase. Results for the 
MIT bag model with m u = mj = 0, m s = 150 MeV and two 
different values of the bag constant B, are presented in the top 
panels whereas those for the Chromo Dielectric (CDM) and 
the NJL models are shown in the bottom ones. The arrow and 
the corresponding number indicates the value (in MeV/fm 3 ) 
of the transition pressure Pq. 
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TABLE I: Stellar properties for the hadronic EoS of Ref. [2?| (Model 1) in the case nuclear matter (only nucleons) or hyperonic 
matter (nucleons+hyperons) and different models for the quark phase. M^ax is the gravitational maximum mass of the pure 
hadronic star sequence, M t hr = M(Po) is the gravitational threshold mass for metastable stellar configurations, M cr (M„) 
the gravitational (baryonic) critical mass, Mfi n is the mass of the hybrid star which is formed by the stellar conversion of the 
hadronic star with M = M cr and assuming baryon number conservation in the process (i.e. Mj in = M^r), finally E conv is the 
total energy, in unit of 10 51 erg, liberated in the stellar conversion. All stellar masses are expressed in units of the solar mass 
M = 1.9889 x 10 33 g. The bag pressure B is given in MeV/fm 
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TABLE II: Same as in Tab. I, but now for the hadronic EoS of Ref. J3]J (Model 2) and three different parametrizations of this 
model: Pi (K x = 236 MeV and x = 1), Pi (K x = 236 MeV and x = 1/3) and P 3 (K^ = 285 MeV and x = 1). 



tion occurs for any of the three parametrizations Pi , Pi 
and P3. In Fig. [2] we plot the Gibbs energy per baryon 
as a function of pressure using: the MIT bag model for 
the Q*-phase and an hadronic parametrization with the 
incompressibility = 236 MeV (Pi) and Koo = 285 
MeV (P3) with x = 1 in both cases. Results for the 
parametrization P2, and the CDM and NJL models are 
not shown for conciseness. It is interesting to note that 
the parameter x, does not influence much the mass and 
radius of the hybrid star maximum mass configurations. 

Note also (see Table [TTJ) that, similarly to Model 1, 
for the parametrizations Pi and Pi , the formation of the 
Q* -phase is possible only in the case of the MIT bag 
model EoS with a low value of the bag constant (B = 85 
MeV fm -3 ). Nevertheless, for the parametrization P3, 
the stable hybrid sequence may be populated from the 
stellar conversion of the critical mass hadronic star if the 
quark phase is described cither with the MIT bag model 
(with B = 85 MeV/fm 3 ) or with the CDM model. As 
already said, no transition is found for the NJL model. 
We also note that in the case of the parametrization P3 
plus the MIT bag model with B = 130 MeV/fm 3 the Q*- 
phase nucleation time at the center of the maximum mass 
(M^ x ) hadronic star is much larger than the age of the 
Universe, and thus, it is extremely unlikely to populate 
the hybrid star branch in this case. In the most favorable 



scenario the possible largest star mass would be 1.624 Mq 
(for CDM). This mass could occur if after, the conversion, 
the star accretes mass from an eventual companion star 
if the object is in a binary system. 



V. SUMMARY AND CONCLUSIONS 

Using the microscopic Brueckner-Hartree-Fock ap- 
proach to describe the EoS of dense hadronic matter we 
have studied the possibility of occurrence of a deconfine- 
ment phase transition into quark matter in neutron star 
cores. Quark matter has been described with three dif- 
ferent models, namely, the MIT bag, the CDM and the 
NJL models. We have concluded that hyperonic hadronic 
stars will not suffer a deconfinement phase transition ex- 
cept if the quark EOS is obtained using the MIT bag 
model with a value of the bag pressure B = 85MeV/fm 3 . 
In this case, however, it is not possible to get a star with 
a mass above 1.54 Mq. On the other hand, we have 
found that if the hadronic matter has no hyperons then 
deconfinement will occur only in very massive stars, with 
M > 2.2Af Q and the stars will decay into a black-hole. 
Within this microscopic approach to the hadronic phase 
the pulsar PSR J1614-2230 would be an hadronic star 
containing only nucleons and leptons. 



5 



We have also studied the possible effect of a hypcronic 
TBF using the model proposed in [3jJ . It was shown that 
for the hardest EOS with = 285 MeV and x = 1 an 
hybrid star could be formed. Only NJL and the MIT 
with B = 130 MeV/fm 3 did not predict a metastable 
star in this case. Within this scenario a maximum mass 
of 1.624 Mq was predicted, very far from the mass of the 
PSR J1614-2230. However, it was also shown that it was 
not so much the hyperonic TBF strcngh but more the 
incompressibility of the nuclconic part of the EoS that 
defines the possible deconfincment transition. 
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FIG. 2: (Color online) Gibbs energy per baryon at zero 
temperature as a function of the pressure for the hadronic 
phase (solid lines) and the Q*-phase (dashed lines). The 
parametrizations Pi (left panels) and P3 (right panels) of the 
Model 2 introduced in [3l( have been used to describe the 
hadronic phase. Results for the MIT bag model with B = 85 
(top panels) and 130 MeV/fm 3 (bottom panels) are presented. 
The arrow and the corresponding number indicates the value 
(in MeV/fm 3 ) of the transition pressure Pq. 
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